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^ ■ ABSTRACT 

As many as 10^ neutron stars populate the Galaxy, but only ~ 10^ are directly observed 
as pulsars or as accreting sources in X-ray binaries. In principle also the accretion of the 

■ interstellar medium may make isolated neutron stars shine, and their weak luminosity could 
, be detected in soft X-rays. Recent ROSAT observations have convincingly shown that neutron 
' stars accreting from the interstellar medium are extremely rare, if observed at all, in contrast 

II . with earlier theoretical predictions. Until now two possible explanations for their elusiveness 

O \ have been proposed: their velocity distribution may peak at ~ 200 — 400 kms~^, as inferred 

^ ■ from pulsar statistics, and this would severely choke accretion; the magnetic field may decay 

^ \ on timescales ~ 10^ — 10^ yr, preventing a large fraction of neutron stars from entering the 

accretor stage. The search for accreting neutron stars has produced up to now a handful of 
promising candidates. While little doubt is left that these objects are indeed isolated neutron 
^ ' stars, the nature of their emission is still controversial. In particular accreting objects can be 

^ , confused with much younger, cooling neutron stars. However, a combination of observations 

and theoretical modeling may help in discriminating between the two classes. 

Subject headings: Stars: Neutron — Stars: X-rays — Stars: Dynamics 



1. INTRODUCTION 

Neutron stars (NSs hereafter) are a well established 
end point of stellar evolution. Although their exis- 
tence was proposed in the thirties, their actual dis- 
covery coincides with that of pulsars and with the 
recognition that pulsars are young, rapidly spinning 
NSs, endowed with a large magnetic field. A few 
years later it was realized that most X-ray binaries, 
and X-ray pulsators in particular, do contain accret- 
ing NSs. 

Altogether observed pulsars and X-ray binaries ac- 
count for ~ 1500 objects. However, as indicated by 
various arguments (see section 4), NSs should repre- 
sent a non-negligible fraction, ~ 1%, of the stars in 



the Galaxy, with a total number as high as 10 . NSs 
are therefore, rather numerous but extremely elusive. 

Let us consider isolated NSs, which supposedly rep- 
resent the bulk of the population. Soon after their 
birth they may show up as radio pulsars, if the beam- 
ing conditions are favorable, and the pulsar phase 
may last some 10^ yr. At formation NSs are hot 
(T ~ 10^^ K), and they may give off thermal radi- 
ation while cooling. Pulsar activity may contribute 
to the warming, but after a few tens of millions years 
the supply of energy, both rotational and internal, 
will be essentially exhausted. From then on, NSs 
are cold, dead objects, apart for some sudden re- 
lease of internal energy, like seismic activity, which 
has been often proposed but never convincingly ob- 
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served. Since the Galaxy is ~ 10^" yr old, one infers 
that the active state of NSs is a tiny fraction of their 
life, or, conversely, that the Galaxy is the graveyard 
of a multitude of dead NSs. 

The dowry of collapsed objects is their gravitational 
field, which in the case of NSs is strong, since their ra- 
dius is about 3 gravitational radii. NSs in the Galac- 
tic halo may become responsible for lensing episodes, 
and in fact they have been invoked in that context 
(see e.g. Wasserman & Salpeter 1994). Also, many 
investigations suggest that the magnetosphere of old 
NSs accreting from dense clouds may be bounded by 
a collisionless shock (see e.g. Arons & Lea 1976), 
although a self-consistent analysis probing its exis- 
tence has not been presented yet. In this picture, 
the Alfven radius could be a possible site for accel- 
eration of cosmic rays (Shemi 1995). Inelastic p — p 
collisions between the accelerated particles and the 
accreting flow may also produce 300 GeV-1 TeV 7- 
rays, that may be detected by Cerenkov experiments 
(Blasi 1996). However, despite these remain promis- 
ing possibilities, up to now no NS has been detected 
unambigously through these techniques. 

Dead neutron stars may, however, be brought to 
life again by accretion from the interstellar medium 
(ISM). This idea, which is the main guideline of 
this review, was proposed some 30 years ago by Os- 
triker, Rees & Silk (1970, hereafter ORS), but for two 
decades it remained latent in the impetuous progress 
of X-ray astronomy. It was only in the early '90s that 
Treves & Colpi (1991) and Blaes & Madau (1993), 
hereafter TC and BM, realized that the ROSAT X- 
ray telescope (launched in May 1990) could have the 
capability of detecting NSs accreting from the ISM. 
These conclusions motivated a renewed interest in the 
field, and several studies aiming to estimate the num- 
ber of detectable sources in the more favorable sites 
in the Galaxy were presented. 

It has become increasingly clear that a realistic es- 
timate of the observability of old, isolated, accret- 
ing neutron stars (ONSs) relies on a thorough under- 
standing of a number of issues, such us 

• the present density and velocity distribution of 
NSs, a subject of interest in itself, strongly related 
to pulsars birth rate and kick distribution. The 
accretion rate is, in fact, inversely proportional to 
the star velocity cubed; 

• the magnetic field and spin evolution of isolated 
NSs, an unsolved problem in compact object as- 
trophysics. Their knowledge is crucial for old neu- 
tron stars, since accretion onto a rotating dipole 
depends critically on the field strength and the ro- 
tation period; 

• the distribution of the interstellar medium (ISM) 



and Giant Molecular Clouds in the Galaxy. The 
ONS luminosity scales linearly with the density of 
the ambient medium, so brighter sources are found 
where the ISM is denser; 

• the spectral properties of radiation emitted by 
neutron stars accreting at low rates, and their de- 
pendence on the magnetic field. 

In all these fields there has been substantial progress 
over the last few years. 

Spectra from accreting NSs for very low luminosi- 
ties have been computed for both unmagnetized and 
magnetized atmospheres by Zampieri et at. (1995) 
and Zane, Turolla & Treves (1999). This regime was 
never explored before, since attention was mainly fo- 
cussed on X-ray binaries, where the luminosity is not 
far from the Eddington limit, millions of times that 
expected from ONSs accreting from the ISM (see sec- 
tion 3). 

The present distribution of NSs has been studied by 
several authors (Paczynski 1990; Blaes & Rajagopal 
1991; BM; Zane et al. 1995) and led to the predic- 
tion that a large number of accreting ONSs should 
be detectable with ROSAT in regions where the ISM 
is denser, such us the closer Giant Molecular Clouds 
(BM; Colpi, Campana & Treves 1993; Zane et al. 
1995), or in the Solar vicinity (Zane et al. 1996b). It 
was also suggested that the collective emission from 
the ONS population could account for part of the 
galactic excess in the soft X-ray background (Zane 
et al. 1995) and for the soft component of the dif- 
fuse X-ray emission from the Galactic Center (Zane, 
Turolla & Treves 1996). 

However, as soon as the first results of the searches 
for isolated NSs were published (Motch et al. 1997; 
Belloni, Zampieri k. Campana 1997; Maoz, Ofek & 
Shemi 1997; Banner 1998a, b), it became apparent 
that theoretical estimates were in excess by a factor 
at least 10-100 over the upper limits implied by ob- 
servations. Still, in the last few years half a dozen 
objects found in ROSAT images appear as good can- 
didates for close-by ONSs accreting from the ISM 
(Stocke et al. 1995; Walter, Wolk & Neuhauser 1996; 
Haberl et al. 1996, 1997; Haberl, Motch k Pietsch 
1998; Schwope et al. 1999; Motch et al. 1999; Haberl, 
Pietsch & Motch 1999). 

The most recent line of research is therefore 
twofold. First it aims to explain the paucity of ob- 
served ONSs. At present, the most likely answers 
to this puzzling question are that NSs have a much 
larger mean velocity than that assumed in the past 
(so the luminosity is drastically reduced), or that the 
magnetic field decays over a timescale ~ 10^ — 10^ 
yr (preventing the star from entering the accretion 
stage). In both cases upper limits on the number of 
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observable ONSs can be used to constrain the popula- 
tion properties (Colpi et al. 1998; Livio, Xu & Frank 
1998; Popov et al. 1999). The second goal is to pro- 
vide a definite assessment of the nature of the seven 
candidates proposed so far. Although their identifi- 
cation with isolated NSs seems rather firm, it is still 
unclear if they are rather young cooling objects, giv- 
ing off thermal radiation at the expense of their in- 
ternal energy, or accreting ONSs. Definite proof may 
come from observations at optical wavelengths, where 
the two classes of sources should show rather distinct 
emission properties (Zane, Turolla & Treves 1999). 

New generation X-ray satellites (CHANDRA, 
XMM, ASTRO E), may either increase the num- 
ber of ONS candidates, or lend further support to 
their paucity. Since these new missions are becoming 
operational, it seems appropriate to review now our 
present knowledge of the subject. 

The structure of the paper is the following. In sec- 
tion ^ we summarize the basic conditions under which 
accretion is possible and section ^ contains a discus- 
sion about the spectral properties. The conditons 
affecting the number of old, accreting NSs are ana- 
lyzed in section ^ Section ^ reviews theoretical in- 
vestigations about the observability of single sources. 
A report of the present status of observations is given 
in section |6| and the observational appearance of ac- 
creting and cooling NSs is discussed in section |^ In 
section ^ possible explanations for the paucity of ac- 
cretors are considered. Our conclusions follow in sec- 
tion ^. 

2. INTERACTION WITH THE INTERSTELLAR 
MEDIUM: VARIOUS TYPES OF NS 

Let us consider a NS of mass M, radius R, spin 
period P, magnetic field B, moving with velocity V 
relative to an ambient medium of number density n. 
There are various possible modes of interaction be- 
tween the NS and the medium (see e.g. Davidson & 
Ostriker 1973; Lipunov 1992, and, specifically for the 
case of isolated NSs, Treves, Colpi &: Lipunov 1993; 
Popov et al. 1999) and these are well illustrated by 
introducing different characteristic lengths. The most 
relevant ones are the accretion, the corotation and the 
Alfven radius, all of them familiar from the theory of 
binary X-ray sources. 

The accretion radius, race defines the region where 
the dynamics of the ISM is dominated by the gravi- 
tational field of the NS and is given by 



2GM 



3 X 10'^ mv^tf cm 



(1) 



where t-io = {V^ + C^)^/^/{W kms"^), m = M/Mq, 
Mp) is the solar mass and C, ~ 10 kms~"^ is the ISM 



sound speed. Equation (|l|) is derived in the frame- 
work of the Hoyle-Bondi theory of accretion. The 



resulting value of race is only approximate, but is suf- 
ficient for our purposes. 

Similarly, the Alfven radius, r^, is the boundary 
inside which the dynamics of the infalling matter is 
dominated by the NS magnetic field. In spherical 
symmetry, 

r A — ~ 

\^/2GMM ) 

~ 2 X lQ^^Bt!,'M-^''R^^"m-^/' cm (2) 



where M is the accretion rate (see eq. |]10[), B12 = 
5/(10i2 G), Mil = M/(10" g s-i) and R^ = R/{m^ 
cm). 

Finally, the corotation radius is obtained by equat- 
ing the angular velocity of the NS with the Keplerian 
angular velocity, 



1/3 



2 X lO^m^/^P^/^ cm . (3) 



In the case of a NS, strong magnetic fields and fast 
rotation may inhibit accretion because of the mo- 
mentum outflow, produced by the spinning dipole, 
and the corotating magnetosphere (ORS; Illarionov 
&: Sunyaev 1975; Davies & Pringle 1981; Blaes et al. 
1992). Three basic conditions must be verified in or- 
der to make accretion possible. First, accretion is 
inhibited if the Alfven radius is larger than the ac- 
cretion radius in which case the system remains in 
the so-called georotator stage. 

Second, at the accretion radius the gravitational 
energy density of the incoming material 



Ug 



GMm„n 



6.5 X 10"^^ Afiir^//^ erg cm"^ (4) 



must be greater than the energy density, Ub-, of the 
relativistic momentum outflow produced by the ro- 
tating i?-fleld. If the fleld is dipolar. 



= [^) (^) (S) ^ 



~ 7.5 X Bi^P R^r^t ergcm""* , (5) 

where ri4 = r/(10"'^^ cm) and rc = cP/2it is the light 
cylinder radius. This condition is met only when the 
NS has spun down to a period 



P ;> Pcrit ~ lQB\!,^M^^'\rA)iARl^^m-'l'' s , (6) 

while during the time in which P < Pcrit the system 
is in the ejector phase, with no accretion occurring. 
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The duration of this phase may largely exceed the 
pulsar life time (~ 10^ yr). 

Since the star is slowing down at the magnetic 
dipole rate, this first barrier is overcome in a typical 

time scale h ~ AB^^^u^"^ Gyr. As noted by BM, 
this value is uncomfortably close to the age of the 
Galaxy. However, if taken literally, it implies that 
a large fraction of ONSs can have spun down suffi- 
ciently, at least if the majority of them are born early 
in the Galactic history. 

After P has increased above Pent, the infalling ma- 
terial proceeds undisturbed until the Alfven radius, 
where the NS magnetic energy density balances the 
matter bulk kinetic energy density. The corotating 
magnetosphere will then prevent the accreting mate- 
rial from going any further, unless the gravitational 
acceleration at the Alfven radius is larger than the 
centrifugal pull, i.e. 



GM 



27r 



rA ■ 



(7) 



This translates into another stronger constraint on 
the value of the period, since it must be 



P;>Pa^ W^B%'M-'/^m-'/^ 



(8) 



(our third condition) otherwise matter will accumu- 
late at the Alfven radius and the system remains in 
the so-called propeller phase. 

The centrifugal barrier at the Alfven radius poses a 
severe problem, since Pa is so large that it cannot be 
reached only by means of magnetic dipole radiation. 
However, there are at least two other effects that may 
play an important role, making accretion possible: 
the decay of the B-field and the torque exerted by 
the accreting material on the NS itself. The accreted 
matter will be spun up by the magnetosphere, and 
it will exert a torque on the NS (see e.g. BM). The 
propeller physics is very complicated, and a detailed 
review of this issue is outside the scope of this pa- 
per. Its thorough understanding requires a full 2- 
D or even 3-D MHD investigation of the interaction 
between accretion flow and rotating magnetosphere 
(see e.g. Toropin et al. 1999). Nevertheless, approx- 
imated expressions for the torque were presented by 
BM (see also, for other spin-down formulae, Lipunov 
& Popov 1995). Based on that treatment, the cor- 
responding spin-down time to Pa turns out to be 



0.04B 



12 



ll/14^-17/28„29/ 



14 



''10 



Gyr, a value adequate to 
allow interstellar accretion even for high values of the 
magnetic field. However, for such high fields accretion 
is probably unsteady and ONSs might appear as tran- 
sient X-ray sources (Treves, Colpi k. Lipunov 1993). 
Numerical simulations indicate that more rapid spin 



down occurs either because of the large mass ex- 
pulsion rate (Toropin et al. 1999; Toropin, private 
communication) or because the material builds up, 
compressing the magnetosphere and becoming unsta- 
ble to large scale mixing with the 5-field (Wang & 
Robertson 1985). 

In any case, if the second barrier is also overcome, 
then matter flows down onto the NS surface and the 
system becomes an accretor. Eventually, at least in 
the simplest picture, if the magnetic field is strong 
enough to channel the accretion fiow, the gas slides 
along the open field lines and the emitting region is 
restricted to two polar caps of radius 

^'^^^ n in3 D-2/7 ,Vl/7 o9/14 1/14 
-Tcav ' ~ ^ ^ 10^12 ^^11 -^6 ^ cm . 

(9) 

The accretion rate is a sensitive function of the star 
velocity and, again within the Bondi-Hoyle theory, is 
given by 



M 



2TT{GMfmpn 

(y2 + C'2)3/2 



lO^^nvfff gs"^ 



(10) 



The total luminosity is then 



L 



GM 



M ~ 2 X 10^^ Mil ergs"^ 



(11) 



which is orders of magnitude below the Eddington 
limit. This implies that matter is very nearly in free- 
fall until it reaches the NS surface. The dynamical 
time is then much shorter than the radiative cooling 
times and virtually no energy is released before the 
fiow hits the outermost stellar layers. Here accreting 
protons are decelerated by Coulomb collisions with 
atmospheric electrons and/or by plasma interactions 
and the fiow stops after penetrating a few 10) 
Thomson depths in the NS atmosphere. The bulk 
kinetic energy of the infalling protons is transformed 
into thermal energy and finally converted into elec- 
tromagnetic radiation. 

An example of the possible stages for a NS 10 Gyr 
old as a function of the star velocity and magnetic 
field is shown in figure |l]. 

3. THE SPECTRUM OF ACCRETING NS 

Due to the combined effects introduced by the re- 
sponse of the detector and absorption by the ISM, 
theoretical estimates on ONS observability and the 
choice of the most favorable energy bands for their 
detection are crucially related to the determination of 
the spectral properties of the emitted radiation and 
to the evaluation of the mean photon energy in par- 
ticular. 
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3.1. Black Body Spectrum 

Many investigators assumed (and still do!) that the 
spectrum emitted by NSs accreting from the ISM is 
a blackbody at the star effective temperature, 



/ L 



1/4 



V47r/i?V 



3.4 X W^lH^Rq 



1/2 



/-i/^K, (12) 



where / is the fraction of the star surface covered by 



accretion and L31 = L/(10 ergs" ) (see eq. 
Although crude, such approximation is basically cor- 
rect, since the density in the deep atmospheric layers 
for T ~ Teff turns out to be 



GMnipTes 



eff 



IQTesmRfiTeff)^^ gcm-3 (13) 



where {Tgff)^ = Tejj/(10^K). For these values of 
density and temperature, the free-free optical depth 
Tff {1^ f f / i^es)Tes IS much larger than unity up 
to hv ^ kT^ff, so thermal equilibrium is estab- 
lished. Compton scattering is not expected to modify 
the spectrum because of the relatively low Thomson 
depth and electron temperature. Cold atmospheric 
electrons emit cyclotron radiation, but for i? ~ 10^ G 
the cyclotron line contribution to the total luminosity 
is vanishingly small and it never exceeds a few per- 
cent even for B ~ 10^^ q (Nelson et al. 1995). These 
considerations show that, at least to first order ap- 
proximation, the emission is peaked at an energy 



E ~ SkTeff 



imLy^'R,"^r'/^ eV 



(14) 



which falls in the soft X-ray range (see ORS). 

For a low magnetic field {B <^ lO^'' G) emis- 
sion/absorption processes (and radiative transfer) are 
much the same as in the non-magnetic case, but the 
emitting area is now substantially reduced because 
<C R. This has an important effect on the emit- 



' cap 



ted spectrum, since, for a given luminosity, the effec- 
tive temperature increases with decreasing / 



T, 



eff 



^eff,B=0 



c-1/4 



47ri?2 



1/4 



~ 13 B^2 M-^^^ Rq m 



1/28 



(15) 



Of course the same argument about the reduced size 
of the radiating region holds also for strongly magne- 
tized NSs {B > lO^'' G), although in the latter case 
the different response properties of the magnetoactive 
plasma must be accounted for. 

The typical values of Tf^jf and L indicate the two 
basic reasons that make the detection of these sources 



extremely difficult: their intrinsic weakness (the lu- 
minosity is orders of magnitude below Solar) , and the 
fact that the their spectrum peaks in the EUV-soft 
X-rays, an energy band which is strongly absorbed 
and not easily accessible even to spaceborne instru- 
mentation. A definite prediction, which is largely in- 
dependent of the details of the spectral distribution, 
is that the X-ray to optical flux ratio must be ex- 
ceedingly large for these sources. 



log 



fx 
fv 



5.5 + 3 log 



eff 



100 eV 



(16) 



for a blackbody spectrum, where fx and fy are the 
fluxes integrated in the intervals 0.5-2.5 keV, and 
3,000-6,000 A, respectively. 

3.2. Beyond the Black Body Spectrum 

The problem of determining the detailed radiation 
spectrum from accreting NSs was first addressed by 
Zel'dovich &; Shakura (1969) and further investigated 
by Alme & Wilson (1973). Both authors considered 
a static, plane-parallel, pure H atmosphere with neg- 
ligible magnetic field, kept hot by the energy released 
by the infalling protons, and focussed on the high ac- 
cretion rates, typical of X-ray binaries. Models are 
characterized by two basic parameters, the accretion 
rate and the penetration depth of protons falling into 
the NS outer layers. Using the same input physics, 
Zampieri et al. (1995) extended the calculations to 
the low luminosity range (L ^ 10^^ erg/s) relevant to 
accretion from the ISM. They found that the emerg- 
ing spectrum is harder with respect to a blackbody 
at Tgff. The hardening factor T^/T^ff, where T-y is 
the radiation temperature, is 2-3 and increases for 
decreasing luminosity. Not surprisingly, the same be- 
haviour is shared by (unmagnetized, H) atmospheres 
around cooling NSs (Romani 1987). In fact, the hard- 
ening is not related to the heat source but to the fre- 
quency dependence of the free-free opacity (the most 
important radiative process). Higher energies decou- 
ple in deeper layers, which are hotter, and the final 
spectrum is, roughly, the superposition of blackbody 
spectra at different temperatures. 

A detailed computation of spectra emerging from 
strongly magnetized, accretion atmospheres has been 
recently presented by Zane, Turolla & Treves (1999). 
Extending previous work, (see e.g. Meszaros 1992 
and references therein) they solved the transfer equa- 
tions in a slab geometry for the two normal modes in 
the (anisotropic) magnetoactive medium coupled to 
the hydrostatic and energy balance equations. They 
have shown that, for L ~ lO^'' — 10^^ ergs~^, the 
hard tail present in non-magnetic models with com- 
parable luminosity is suppressed and the X-ray spec- 
trum, although still harder than a blackbody at T^f / , 
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is nearly planckian in shape. At X-ray energies the 
spectral distribution closely resembles that of a cool- 
ing, hydrogen atmosphere with similar properties (see 
e.g. Shibanov et al. 1992 and section |^. However, 
the two models (accreting/cooling) differ rather sub- 
stantially at optical wavelengths, where the former 
exhibits a definite excess over the Rayleigh-Jeans tail 
of the X-ray best-fitting blackbody (see figure §). A 
more thorough comparison of the properties of ac- 
creting and cooling models, and of the observational 
relevance of the optical excess as a signature of ac- 
cretion, is postponed to section |^. 

In concluding this section, we remark that, al- 
though this has little effect on the bulk of the emis- 
sion, for B ~ 10^^ G the (broad) electron cyclotron 
line, centered at ~ 11.6 i?i2 keV, becomes a promi- 
nent spectral feature over the falling-off Wien con- 
tinuum, as pointed out by Nelson et al. (1995). The 
detection of the cyclotron line in conjunction with 
polarization measures near the proton cyclotron en- 
ergy (where the degree of polarization crosses zero) 
may prove a powerful tool in determining the field 
strength, even in the absence of pulsations (Zane, 
Turolla & Treves 1999; Pavlov & Zavlin 1999). 

4. CONDITIONS AFFECTING THE NUMBER OF 
ACCRETING ONS IN THE GALAXY 

4.1. The Number of NSs in the Galaxy 

The total number of neutron stars in the Galaxy 
is still uncertain, varying between Ntot ~ 10® and 
~ 10^. These figures are inferred from constraints on 
the nucleosynthetic yields produced in core collapse, 
and from estimates on the present rate of Type II su- 
pernova events, as observed in nearby spiral galaxies. 

According to Arnett, Schramm & Truran (1989), 
core collapse events occur at a rate of one every 10 
years and, when averaged over the entire Galaxy life- 
time, can account for nearly half of the Galactic Fe 
abundance. This would imply ~ 10^ NSs if the disk of 
the Milky Way has a mass of ^ W^^Mq and an age 
of ~ lO^'^ yr. Since thermonuclear events involving 
the explosion of C/0 white dwarfs may contribute to 
the Fe yield, Ntot cannot be determined with higher 
accuracy. 

The present rate of Type II supernova explosions, 
derived from Capellaro et al. (1997; see also Madau, 
Delia Valle & Panagia 1999) is ~ 0.5 HIq SNu, where 
1 SNu is 1 supernova per 100 yr per IO^^Lbq and 
/i50 the Hubble constant in units of 50 kms~^Mpc~^. 
For a Galactic blue luminosity of 5 x IO^^^L^q, this 
estimate implies a supernova event every 30 years, 
approximately. Thus, Ntot ~ 3 x 10® NSs should be 
present in the Milky Way. However, this estimate 
may just give a lower limit for Ntot since the Galaxy 
may have experienced phases of enhanced star forma- 



tion at early epochs. The rate in the past could have 
been higher by a factor of 10 (Madau et al. 1996) 
yielding again Ntot close to 10^. 

4.2. Space and Velocity Distribution of NSs 

The visibility of accreting ONSs is severely con- 
strained by the present velocity distribution, in par- 
ticular by the extent of a low velocity tail. As old NSs 
outnumber young objects, in all likelihood active as 
pulsars, potential accretors belong predominantly to 
the dynamically old population which had the time 
to evolve in the Galactic potential. The velocity dis- 
tribution of ONSs and their spatial density in the 
volume accessible to observations can be inferred sim- 
ulating a large, statistically significant number of or- 
bits in the Galaxy. The strategy followed by a number 
of authors (BM; Blaes & Rajagopal 1991; Frei, Huang 
& Paczyiiski 1992; Zane et al. 1995) is to start with 
a seed population, born in the disk with natal kicks 
determined according to a specific velocity distribu- 
tion, that evolves dynamically in the Galactic poten- 
tial. Customarily the distribution of the kick velocity 
is approximated as isotropic and Gaussian. The neu- 
tron star birth rate and spatial distribution mirror 
those of massive stars in the Galaxy, while pulsars, 
representative of the parent population, are used to 
constrain the natal velocity dispersion. 

The first to use available velocity information in 
the statistical analysis of young NSs were Narayan &: 
Ostriker (1990). They found that observations of pe- 
riods and magnetic fields from a sample of about 300 
pulsars are well fitted by invoking the presence of two 
different populations of NSs at birth, slow (S) and fast 
(F) rotators. Both have an isotropic Gaussian distri- 
bution (relative to the circular speed), the F rota- 
tors being characterized by a mean velocity ~ 60 
kms~^, and by a scale height zp ~ 150 pc, much 
lower than the S ones for which Vs ~ 240 kms~^ 
and zs = 450 pc. From this fit, the evolved distri- 
bution function of ONSs has been obtained by Blaes 
& Rajagopal (1991), BM, and Zane et al. (1995). 
As far as accretion onto ONSs is concerned, only the 
F population, which represents ~ 55 % of the total 
number, is relevant. Results from Zane et al. (1995) 
show that, after secular evolution, in the local region 
7.5 kpc < R < 9.5 kpc, ONSs are characterized by 
a mean scaleheight {z) ~ 250 pc and by a mean ve- 
locity, averaged over \z\ < 200 pc, {V) ~ 78 kms~"^. 
The number density of NSs, within 2 kpc from the 
Sun, turns out to be no ~ 3 x 10'"^ {Ntot /W^) pc~^. 

The previous scenario was challenged by Lyne & 
Lorimer (1994), who suggested that NSs acquire at 
birth velocities significantly higher than those of both 
the F and S sub-populations of Narayan & Ostriker, 
the imprint of an anisotropic core collapse during 
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a type II supernova. These authors derived from 
the observations a typical mean velocity of about 
450 kms~^. Recently, on completion of a quantita- 
tive analysis, the kinematics of radio pulsars was as- 
sessed using population synthesis models and a sam- 
ple of pulsars with improved statistics and revised 
distance scale. Yet, the results depicting the true 
pulsar population are still not univocal, due to dif- 
ferent criteria used in selecting the sample and in 
treating observational errors. While Lorimer, Bailes 
8z Harrison (1997) find a mean natal velocity of 
500 kms^^, Hansen Sz Phinney (1997) claim a lower 
value, ~ 250-300 kms"^ Cordes & Chernoff (1998) 
suggest a possible decomposition of the velocity dis- 
tribution into two components with mean velocities of 
197, 700 kms~"^ and fractions 0.84, 0.16, respectively. 
Despite these controversies, recent observations seem 
to point towards an increase in the mean velocity of 
young NSs relative to Narayan & Ostriker. 

Accreting NSs belong necessarily to the low velocity 
tail of the population, and at present the low statis- 
tics prevents a determination of the extension of such 
a tail with the required accuracy. Hartmann (1997), 
Hartmann et al. (1997) and Iben & Tutukov (1998) 
do not exclude the presence of a prominent low veloc- 
ity component. Dynamical heating (Madau & Blacs 
1994) may influence the evolution of slow NSs. This 
process, observed in the local disk stellar population, 
causes the velocity dispersion to increase with time as 
a consequence of scattering by molecular clouds and 
spiral arms (Wielen 1977). If ONSs participate in the 
same process, dynamical heating over the lifetime of 
the Galaxy may scatter a fraction of low velocity stars 
to higher speeds and this, in turn, could devoid the 
distribution of slow stars. 



4.3. The Distribution of the Interstellar Medium 

The discussion of the properties of accreting 
sources, and the study of the detectability of ONSs, 
require a detailed mapping of the interstellar medium. 
For a given source, the interstellar density determines 
both the intrinsic luminosity and the amount of ab- 
sorption along the line of sight. The situation is more 
uncertain and complex when dealing with particular 
regions such as the center of the Galaxy or the So- 
lar proximity, while the gas structure is smoother on 
larger scales and numerical fits to the gas distribution 
have been published in the literature (see e.g. Dickey 
&: Lockman 1990 for a detailed review). 

In particular, observational data from Lya and 21 
cm absorption measures show that the ISM distribu- 
tion of both cold and warm HI is nearly constant in 



radius while its z-dependence can be fitted by 

riHi = ^^n^exp +n3exp (^-j^ (17) 

with m = 0.395, n2 = 0.107, = 0.064; ai = 212, 
(T2 = 530 and h = 403 pc; Ui and ai are in cm~^ and 
pc, respectively. The applicability of the previous ex- 
pression is restricted to the range 0.4 < R/Rq < 1, 
where R is the galactoccntric radius and Rq = 8.5 kpc 
is the distance of the Sun from the Galactic center. 
The gas layer has a scale height of about 230 pc in the 
vicinity of the Sun, while for R < O.^Rq it shrinks to 

100 pc and in the outer Galaxy it expands linearly 
up to ~ 2 kpc. The other important contribution to 
the total ISM density comes from molecular hydro- 
gen. The best tracer of H2 is the CO molecule, and 
observational data suggest a local gaussian distribu- 
tion with a scale height of ~ 60 — 75 pc. Observa- 
tions, however, are much less conclusive as far as the 
midplane density is concerned (Bloemen 1987), also 
because it may significantly depend on R (Dc Boer 
1991). As a first approximation, the H2 distribution 
can be modeled as a gaussian with central density 0.6 
cm~3 and FWHM 70 pc (De Boer 1991). The ion- 
ized component gives only a very small contribution 
to the total density, and can be neglected. 

The evidence for a complex sub-parsec structure 
in the ISM has been recently found by a number of 
authors on different scales, from 10^ — 10^ AU down 
to 10 - 100 AU (see e.g. Meyer & Lauroesh 1999 
and references therein). On all the sampled scales 
observations imply dense concentrations of atomic 
gas {uh ^ lO"' cm~^) in otherwise diffuse sight fines. 
These local condensations can not be accounted for 
in the standard McKee &: Ostriker (1977) equilibrium 
pressure model for the ISM, thus further interstellar 
line mapping is required to increase the small scale 
sky coverage and to probe the sophisticated spatial 
structure. On the other hand, their existence could 
explain why some ONS candidates are found in very 
low average density regions. 

Even leaving aside the issue of the ISM dumpiness 
on small scales, the structure of the ISM at distances 
^ 1 kpc the ISM is highly anisotropic and its com- 
plex morphology and physical state are currently a 
very active field of research. At the same time, this 
represents a particularly favorable site, containing a 
high number of ONSs that, due to their vicinity, may 
be relatively easy to detect. Assuming a spatial den- 
sity of ~ 3 X 10""^ pc"-'^ (section 4.2, BM, Zane et 
al. 1995), about 150 ONSs are in fact expected in a 
sphere of radius 50 pc centered on the Sun. Quite dis- 
appointingly, the observability of these close-by NSs 
as accretion powered sources is severely hindered by 
the shortage of fuel. In fact, the Sun is surrounded 
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by a region, the Local Bubble, where the plasma 
has both very low density {n ~ 0.05 — 0.07 cm~'^) 
and high temperature (T ^ 10^ K). In the scenario 
proposed by McKee & Ostriker (1977), see also Cox 
& Anderson (1982), Cox (1983), the hot gas would 
fill ~ 70 — 80% of the interstellar space and a large 
number (~ 2 x lO'^) of cool (T ~ 80 K), roughly 
spherical clouds are expected to be present. Obser- 
vational data support this model for the region be- 
yond ~ 50 — 100 pc from the Sun (Knude 1979), but, 
as discussed by Paresce (1984), soft X-ray, radio and 
color excess surveys seem to indicate that no clouds 
are present at smaller distances and that the denser 
material is more probably organized into large, elon- 
gated, moving fronts within ~ 50 pc. The Sun it- 
self is embedded in a medium, the Local Fluff, which 
is warm (T ~ 10'^ — 10^ K) and slightly overdense 
(n ~ 0.1 cm~^) with respect to the Local Bubble on 
scales ^ 20 pc (Diamond, Jewell & Ponman 1995). 

The present picture indicates that the contour of 
the Local Bubble in the Galactic plane is highly asym- 
metric, with four major discontinuities in four differ- 
ent Galactic sectors (Paresce 1984). In particular, a 
wall of neutral hydrogen is located very close to the 
Sun in the second quadrant, 15° < I < 120°. Accord- 
ing to Paresce (1984), the wall is roughly parallel to 
the I = 330° — 150° axis and is located at d < 16 pc, 
with an estimated depth of about 35 pc. The Nhi 
contours presented by Frisch &: York (1983) are gen- 
erally farther away, with the denser material (n ~ 1 
cm~^) at ~ 90 pc from the Sun. Welsh et al. (1994) 
have derived a highly asymmetric contour of the Lo- 
cal Bubble that in the second quadrant is roughly in- 
termediate between those presented by Paresce and 
Frisch & York. The minimum radius of the local cav- 
ity has been estimated to be ~ 25 — 30 pc, but, as 
stressed by the same authors, their indirect method 
could produce an underestimate of Nhj at distances 
smaller than 50 pc. An analysis of ROSAT EUV data 
(Diamond, Jewell & Ponman 1995) has shown that n 
reaches ~ 1 cm~^ at ~ 25 — 30 pc and this result 
seems to be in agreement with the asymmetric con- 
tour found by Welsh et al. more than with those of 
Paresce, Frisch &: York (see also Pounds et al. 1993). 

5. ARE ONS DETECTABLE ? 

As we have seen in the previous sections, the ex- 
pected luminosity from an ONS accreting from the 
ISM is ^ 10^^ ergs^i with a mean photon energy in 
the interval ~ 0.1 — 1 keV. The expected flux at Earth 
is <^ 10~^^(d/100 pc)~^ ergcm~^s~^ (where d is the 
distance) which translates into ^ 1 ROSAT PSPC 
counts"^, well above the typical threshold of the 
ROSAT All Sky Survey (RASS, ~ 0.02 count s'^). 
The sensitivity of the ROSAT PSPC, together with 



its unprecedented capability in the soft X-ray band 
(~ 0.1 — 2 keV) put accreting ONSs within obser- 
vational reach. At the same time the ROSAT HRI 
had a very high angular resolution, which is decisive 
for follow-up pointed multiwavelength observations 
of detected sources. 

5.1. Selection Criteria for isolated NSs 

Since accreting ONSs should be detectable with 
ROSAT, it is important to establish the most rele- 
vant selection criteria. The basic ones, mostly in- 
dependent of details about the spectral shape, were 
already proposed by TC: 

• the X-ray to optical flux ratio should be extreme 



(see section 3T), and the source should show no 
emission in other energy bands (radio, IR, 7-rays); 

• the X-ray spectrum should be thermal and soft, 
with typical temperatures T « 100 eV (see sec- 
tion p]) ; 

• the source must be weak, L ^ 10'^^ erg/s; 

• the low luminosity makes only close-by sources 
{d PS 100 pc) potentially detectable. This, in 
turn, implies that the column density Njj must 
be ^ 10^^ cm~2; 

• since M oc n, the spatial distribution of sources 
should correlate with the denser regions of the 
ISM; 

For close-by sources the visual magnitude should 
be ~ 25-26 and this might make both the proper 
motion and the parallax measurable, especially if iso- 
lated NSs are fast objects. Would the velocity dis- 
tribution be rich in slow stars, more distant sources, 
that may appear in Deep Exposures, must concen- 
trate along the Galactic plane because the scale 
height of both gas and ONSs is ~ 200-300 pc. 

5.2. Number of Detectable Sources with ROSAT 

The observed count rate for a source of monochro- 
matic luminosity L^, is given by 



CR 



CO (i 1/ 

exp -(a^y/V//)^!^!— counts"-^ (18) 
47ra^ nv 



where ay is the interstellar absorption cross section 
and Ajj is the detector effective area. \i CR repre- 
sents the detector sensitivity limit, the previous ex- 
pression can be used to find the maximum distance, 
dmax , at which a star of given luminosity L produces 
a count rate above the threshold. The total number 
of detectable ONSs in a given field is obtained by 
integrating the NS distribution function, f{r,l,b,v), 
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over the solid angle subtanding the field and out to 

dmax 

/f fdmax 
dv dn f{r,l,b,v)dr, (19) 

Jn Jo 

where I, b are the usual Galactic coordinates. 

Using the distribution function of Paczynski 
(1990), TC concluded that ~ GOOOiVg sources, lo- 
cated within 500 pc, should appear in the RASS. 
They assumed polar cap accretion {B = 10^ G), a 
blackbody spectrum and considered a low value for 
the uniform ISM density (n = 0.07 cm~'^), appropri- 
ate to the local bubble. A more thorough investiga- 
tion of ONS detectability has been presented by BM. 
They followed the evolution of the population by in- 
tegrating the orbits of 10^ stars moving in the Galac- 
tic potential over the Galaxy lifetime. Newly born 
NSs were assumed to belong to the F sub-population 
of Narayan & Ostriker (1990). In the solar proxim- 
ity the present ONS mean velocity was found to be 
~ 80 kms^^ and the fraction of stars with y ^ 40 
kms^^~ 20%, two times lower than what is assumed 
by TG. Their results were, nevertheless in substan- 
tial agreement with those of TC and predicted that 
several thousands sources should be present in PSPC 
exposures. Results by BM were improved by Man- 
ning, Jeffries & Willmore (1996), who computed the 
number of possible detections in the Rosat Wide Field 
Gamera survey with a Monte Garlo simulation. By 
taking into account effects of a varying magnetic field, 
NS mass, and exposure times during the survey, they 
were able to predict a significant depletion in the pre- 
dicted numbers of detections if the magnetic field has 
not decayed to values <^ 10^^ G. 

Relatively high count rates are expected from ONSs 
accreting in the denser regions closer to the Sun. Of 
course, these regions must be large enough to con- 
tain a statistically significant number of slow NSs 
which are the most luminous. In this respect, dense 
molecular clouds, where the ambient density is ~ 100 
times higher than the average, provide a very favor- 
able environment for observing ONSs. Estimates by 
BM, Golpi, Gampana & Treves (1993) and Zane et al. 
(1995), using different assumptions about the emit- 
ted spectrum, indicated that the closer clouds should 
harbor some relatively luminous sources. A detailed 
analysis of the observability of ONSs accreting in the 
closest overdense regions of the solar neighbourhood 
has been presented by Zane et al. (1996a). As they 
pointed out, although the local interstellar medium is 
underdense and relatively hot, it contains at least one 
region, the Wall, in which the gas density is ~ 1 cm~^. 
The Wall extends between ~ 20 and 50 pc and its 
angular size is ~ 6000 deg^. According to these au- 
thors, about 10 ONSs are expected to be detectable 



in the Wall directions and, due to their vicinity, these 
sources should appear at the relatively high flux limit 
of 0.1 counts/s in the PSPG survey. 

Besides the detectability of individual sources, 
ONSs could also reveal themselves through their con- 
tribution to the diffuse X-ray background. This 
point, originally suggested by ORS, was addressed by 
BM, who found that the integrated flux above 0.5 keV 
is negligible under the hypothesis of blackbody emis- 
sion. Zane et al. (1995) have readdressed this issue, 
considering polar cap emission and using synthetic 
spectra, and concluded that magnetized, accreting 
ONSs can contribute up to 20% of the unresolved 
soft X-ray excess observed at high latitudes in the X- 
ray background (see also De Paolis Sz Ingrosso 1997 
for a discussion of the possible contribution from NSs 
in the Galactic halo). 

Since the diffuse emissivity of accreting ONSs de- 
pends on their spatial concentration and on the den- 
sity of the ISM, the Galactic Genter, where both 
the star and gas density are very high, could pro- 
vide an excellent site for revealing ONS emission. 
The GG is a well known source of diffuse emission, 
investigated by many missions and in particular by 
GRAN AT. ART-P data in the 2.5-30 keV band 
(Sunyacv, Markevitch & Pavlinsky 1993; Markevitch, 
Sunyaev & Pavlinsky 1993) show the presence of an 
elliptical diffuse source with different properties be- 
low and above ~ 8 kcV. Assuming that ONSs repre- 
sent ~ 1% of stars in the GG and that their distribu- 
tion in both physical and velocity space follows that 
of low-mass stars, Zane, Turolla & Treves (1996) were 
able to reproduce satisfactorily the GRANAT data at 
the lower energies. 

6. ROSAT RESULTS: THE CRUDE REALITY OF 
OBSERVATIONS 

The search for old, isolated NSs in ROSAT images 
is a by-product of "complete" surveys. Typically one 
defines a flux threshold and lists all sources which 
arc brighter than that. The threshold should be such 
that the probability of confusion with the background 
noise is small. The second step is the search for an 
optical counterpart on the basis of positional coin- 
cidence and peculiar colours. Possible optical can- 
didates are then studied spectroscopically to see if 
the source belongs to one of the known classes of X- 
ray emitters. All sources for which no viable optical 
counterpart was found are classified as NOIDs (non- 
optically identified objects). Isolated NSs candidates, 
for which no counterpart is expected down to a visual 
magnitude ~ 26, must be looked for among NOIDs. 

Generally NOIDs are further studied by reducing 
their X-ray error box. The lack of obvious coun- 
terparts in archival optical images, like the POSS, 
is used to set a lower limit on the X-ray to optical 



10 



Isolated Neutron Stars 



flux ratio. Sources with fx/ fv ^ 1000 are classified 
as potential isolated NS candidates. Obviously, the 
number of NOIDs in a given region is an absolute 
upper limit on the number of ONSs. 

RASS and deep field exposures of the most favor- 
able sites (such as GMCs) have been searched for ac- 
creting ONSs and, although the various surveys were 
conducted with rather inhomogeneous criteria, an im- 
portant result soon became apparent. The upper lim- 
its on ONS candidates were substantially below ex- 
pectations estimated by, at least, one order of magni- 
tude. Results from the surveys are reviewed in section 
6.5 and the implications of this (negative) conclusion 
are examined in section ^. 

Despite ONSs being much more elusive than orig- 
inally expected, the observational efforts of the last 
few years have produced some promising candidates 
and their number keeps increasing: one in 1995, two 
in 1996, one in 1998 and three in 1999. All of them 
were discovered, or re-discovered, as in the case of 
MS 0317.7-6647 and RXJ 185635-3754, with ROSAT, 
both serendipitously and in survey data (mainly the 
ROSAT Bright Survey, RBS). The main properties of 
the seven candidates discovered so far are reviewed 
in the following subsections and summarized in ta- 
ble |. 

Although present observations leave little doubt 
that these sources are associated with isolated NSs, 
the exact nature of the mechanism powering their 
emission is still controversial: either accretion of the 
ISM onto old neutron stars or thermal radiation from 
middle-aged, cooling NSs. As we discuss in section 0, 
the emission properties of these two classes of objects 
are rather similar and the present data cannot pro- 
vide convincing evidence in favor of either possibility. 

6.1. MS 0317.7-6647 

This source was the first to be proposed as an ac- 
creting ONS, and, although its association with a 
neutron star was never supported by further obser- 
vations both in X-rays and in the optical, we include 
it in our list. 

MS 0317.7-6647 was already present in the Ein- 
stein Medium Sensitivity Survey ( Gioia et al. 1990) 
and then re-observed with ROSAT by Stocke et al. 
(1995). It is a soft, very weak (~ 0.03 count s~^) 
source in the field of the nearby galaxy NGC 1313. 
Prom the three available X-ray images there is evi- 
dence that this source is variable over the time scale 
of years. The only possible optical counterpart has 
my = 20.8, which brings the X-ray to optical flux 
ratio to fx/fv ^ 60. The simultaneous lack of radio 
emission seems to rule out the possibility of an ex- 
treme BL Lac object leaving only a compact object 
as a reasonable option. Identifying the source with 



an X-ray binary in NGC 1313 leads to a luminosity 
L ~ lO^'^ ergs~^ and requires a mass of the compact 
object ~ 5OM0. Such a large value for the black hole 
mass, and the high Galactic latitude of the source, 
which makes it unlikely to be an X-ray binary in 
our own galaxy, were taken as clues in favor of the 
isolated NS option. The PSPC spectrum, which is 
well fitted by a blackbody with T ~ 200 eV, and 
the derived luminosity, L ~ 1.7 x 10^'^ erg/s placing 
the object at 100 pc, are fully compatible with those 
of an isolated, nearby neutron star. The presence 
of a galactic IR cirrus cloud in the field of view of 
the source supports this possibility, while the sign of 
variability argues against the hypothesis of a cooling 
neutron star (Brazier & Johnston 1999). 

6.2. RXJ 185635-3754 

RXJ 185635-3754 is a rather bright source, already 
present in the EINSTEIN slew survey, and observed 
with ROSAT by Walter, Wolk k Neuhauser (1996). 
The PSPC count rate, ~ 3.6 counts"^, corresponds 
to an X-ray luminosity L ~ 5 x 10^^ erg s~^ for a dis- 
tance of 100 pc. The observed spectrum is very soft 
and the best fit with a blackbody gives a temperature 
of about 60 eV. The original claim that this source 
is an isolated neutron star was supported by the ab- 
sence of an optical counterpart down to my ~ 23, 
giving fx/fv > 7000. 

Further X-ray and optical observations were per- 
formed by Neuhiiser et al. (1997) and Campana, 
Mereghetti & Sidoli (1997) in an attempt to iden- 
tify the optical counterpart and to assess its distance. 
The source appears projected against the molecular 
cloud R CrA and the column density derived from 
X-ray data {Nh ~ 2 x 10^° cm~^) is less than that in 
the cloud itself. The cloud is estimated to be ~ 120 
pc away, so the source is likely to be closer. X-ray 
spectral fits with more sophisticated models of cool- 
ing and accreting NS atmospheres were presented by 
Pavlov et al. (1996) and Campana, Mereghetti & 
Sidoli (1997). Unmagnetized cooling spectra with dif- 
ferent chemical compositions predict an optical flux 
that is either too large or too small, while both mag- 
netized cooling and unmagnetized accreting hydrogen 
models are consistent with the data only if the source 
is as close as ~ 10 pc, and this seems difficult to rec- 
oncile with the derived values of the column density. 

The real nature of RXJ 185635-3754 became even 
more puzzling when Walter & Matthews (1997) re- 
ported the discovery of a possible optical counterpart 
with HST. Their candidate, selected on the basis of 
positional coincidence and its blue colour, is a faint, 
stellar-like object with my ~ 25.6, giving an X-ray 
to optical fiux ratio ~ 75000. They determined the 
optical fiux at two wavelengths, A = 3000,6060 A, 



Treves et al. 



11 



Table 1 
Isolated NS candidates 



Source 


PSPC 


Tbb 


Nh 


log/x//y 


Period 


Refs. 




count 


eV 


1020 cm-2 




s 




MS 0317 


0.03 


200 


40 


> 1.8 




a 


RX J 1856 


3.64 


57 


2 


4.9 




b, c, d, 


RX J0720 


1.69 


79 


1.3 


5.3 


8.37 




RBS1223 


0.29 


118 


~ 1 


> 4.1 




i 


RBS1556 


0.88 


100 


< 1 


> 3.5 




i, 1 


RX J0806 


0.38 


78 


2.5 


> 3.4 




m 


RX J0420 


0.11 


57 


1.7 


> 3.3 


22.7 


n 



Note.— a=Stocke et al. (1995); b=Walter, Wolk & Neuhauser (1996); c= 
Neuhauser et al. (1997); d=Campana, Mereghetti &; Sidoli (1997); e= Walter 
& Matthews (1997); f=Haberl et al. (1997); g=Motch & Haberl (1998); h = 
Kulkarni & van Kerkwijk (1998); i= Schwope et al. (1999); 1 = Motch et al. 
(1999); m = Haberl, Motch & Pietsch (1998); n = Haberl, Pietsch & Motch 
(1999) 



and both points lie above the Rayleigh- Jeans part of 
the blackbody spectrum which best-fits the X-rays. 
The presence of a definite optical excess, about a fac- 
tor 3, points toward a more complex spectrum and 
no definite explanation for this has been found up to 
now (see section ^ for details). 

6.3. RXJ 0720.4-3125 and RX J0420.0-5022 

RXJ 0720.4-3125 was observed for the first time 
with ROSAT by Haberl et al. (1996, 1997). The 
source was detected with the PSPC at a level of ~ 1.7 
counts"^ and the soft spectrum is well fitted by a 
blackbody at T ~ 80 eV. The column density de- 
rived from the fit is low, Nh ~ lO^^ cm"^^ placing the 
source at a relatively close distance, d ^ 100 pc, with 
a luminosity L ~ 2.6 x 10^^((i/100pc)2 ergs~^. Pre- 
liminary optical searches produced a limiting magni- 
tude for the counterpart of mv > 21, implying a flux 
ratio > 500. 

Although the X-ray properties make RXJ 0720.4- 
3125 and RXJ 185635-3754 quite similar, RXJ 
0720.4-3125 has been found to pulsate in X-rays, with 
a period P = 8.39 s. If indeed this source is pow- 
ered by accretion, the knowledge of the luminosity 
and of the period allows an estimate of the magnetic 
field strength, as already pointed out by Haberl et 
al. (1997). The condition that the corotation radius 
is larger than the Alfven radius (see eq. Q) implies 



B ^ 10 G, indicating a substantial decay of the 
magnetic field in this object. 

Haberl et al. (1997) suggested that such values 
of B and P may be the result of common enve- 
lope evolution in a high mass X-ray binary. In this 
case the history of RXJ 0720.4-3125 would bear re- 
semblance to that of the Anomalous X-ray Pulsars 
(AXPs; Mereghetti & Stella 1995) if the last are inter- 
preted as accreting from a residual circumstellar disk 
(van Paradijs, Taam &: van den Heuvel 1995). Their 
observed periods, in the 5-10 s range, are surpris- 
ingly similar to that of RXJ 0720.4-3125. Konenkov 
& Popov (1997) and Wang (1997; see also Colpi et al. 
1998) have shown that the present position of RXJ 
0720.4-3125 in the B — P plane is consistent with the 
long term evolution (~ 10^ — 10^^ yj-'^ gf a slow (~ 20 
kms~^) NS, born with canonical values of the mag- 
netic field and period, which experienced field decay 
over a timescale ^ 10^ yr. The two evolutionary sce- 
narios (common envelope or spontaneous decay) will 
lead to the same kind of object, i.e. a slow, accret- 
ing NS, but on quite different timescales, ^ lO"^ vs. 
109 - 10^0 yr. 

An alternative explanation to the observed prop- 
erties of RXJ 0720.4-3125 was recently proposed by 
Heyl & Kulkarni (1998; see also Wang 1997 and Heyl 
&: Hernquist 1998). In this picture the source is pow- 
ered by the release of internal energy, a magnetar 
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kept hot by magnetic field decay. For B ~ 10 G 
the time required to cool down the NS to the ob- 
served temperature is ~ 10^ yr, so contrary to the 
accretion scenario, RXJ 0720.4-3125 should contain 
a young NS with a decaying field. Heyl Sz Kulkarni 
argued that the source could be the descendant of 
an anomalous X-ray pulsar. In this interpretation, 
AXPs are a class of isolated neutron stars with un- 
usually high magnetic fields, powered by its decay 
(Thompson & Duncan 1996; see also Colpi, Geppert, 
& Page 1999). 

Recent follow-up optical observations (Motch & 
Haberl 1998; Kulkarni & van Kerkwijk 1998) led to 
the identification of a possible optical counterpart, 
a faint, bluish object with ms ~ 26.6. Again the 
identification is based on positional coincidence and 
colour and gives fx/fv ~ 2 x 10^. The observed 
optical flux, similarly to RXJ 185635-3754, exceeds 
that predicted by the blackbody spectrum which fits 
the ROSAT data by a factor ~ 5 (Kulkarni & van 
Kerkwijk 1998). 

RX J0420.0-5022 is the last isolated NS candidate 
discovered so far and the only one, together with RX 
J0720, to exhibit a modulation in the X-ray flux with 
a period ~ 22.7 s (Haberl, Pietsch & Motch 1999). 
It has the lowest observed count rate (~ 0.11 PSPC 
count s~^), but shows all the distinctive features com- 
mon to other sources in this class: a thermal spec- 
trum with T ~ 57 eV, iV/f ~ 2 X 10^° cm'^ which 
give L ~ 2.7 X 1030(d/100pc)2 ergs'^ NTT optical 
observations revealed no unusual object in the X-ray 
error box down to a magnitude m,B ~ 25, impliying 
an X-ray to optical flux ratio larger than ~ 20000. 

6.4. RBS1223, RBS1556 (RX J1605.3+3249) and 
RX J0806.4-4132 

RBS1223 and RBS1556 are two rather bright 
sources, ~ 0.3 and ~ 0.9 PSPC counts"^ respec- 
tively, detected in the ROSAT Bright Survey (RBS) 
and then observed in follow-up pointings by Schwope 
et al. (1999) and Motch et al. (1999). Both 
sources have similar properties, a thermal X-ray 
spectrum with little absorption (T ~ 100 cV, Nh ~ 

1020 

cm ). Optical images taken at Keck show no 
possible optical counterpart down to tub ~ 26 for 
RBS 1223 (Schwope et al. 1999) and to my 24.2 
for RBS1556 (Motch et al. 1999), giving a lower limit 
for the X-ray to optical flux ratio in excess of 10000. 

RX J0806.4-4132 has been serendipitously discov- 
ered by Haberl, Motch & Pietsch (1998) in archival 
PSPC observations of a candidate supernova rem- 
nant. The quite large angular distance from the 
source makes a physical association unlikely. The 
source has all the distinctive features of the other 
candidates, the spectrum is thermal with T ~ 80 eV 



and the column density is low, Nh ~ 2.5 x 10^^ cm~^. 
Preliminary optical searches revealed no unusual ob- 
ject in the X-ray error box down to a blue magnitude 
~ 24. 

6.5. Results from ROSAT Surveys 

Belloni, Zampieri &; Campana (1997) have per- 
formed a systematic search for ONSs in two molecu- 
lar clouds in Cygnus, the Rift and 0B7, by analyz- 
ing archive ROSAT PSPC pointings; this is the only 
deep field survey available up to now. Observations 
cover only a small fraction of the clouds, ~ 5%, and 
contain 109 sources above a flux limit of ~ 0.0015 
counts"^. For 105 of them an optical counterpart 
was identified leaving 4 NOIDs, with no counterpart 
above ma ~ 20. Although presently inferred values 

fx/ fv ~ 1 are not large enough to qualify these 
sources as strong ONS candidates, their positional co- 
incidence with dense clouds together with their hard- 
ness ratio ^ 0.3 make them worthy of future inves- 
tigations. Belloni, Zampieri & Campana estimated 
the number of detectable ONSs in the searched area 
following the criteria of section 5.1, and found that 
it exceeds the number of NOIDs by a factor 3-10. 
Since ONSs can account only for a fraction of NOIDs 
this results provides a clear indication that theoreti- 
cal figures produced so far overestimated the actual 
number of detectable sources. A similar conclusion 
follows from the comparison between the predicted 
number of resolved sources, ~ 10 deg~^ and that of 
NOIDs in high-latitude PSPC pointings, 30 deg"^ 
discussed by Zane et al. (1995) in connection with 
the possible contribution of ONSs to the diffuse X- 
ray background. As they noted, this would imply 
that one NOID in three should be an ONS. 

A more complete analysis of a larger sample region 
in Cygnus, ~ 64.5 deg^, comprising a large part of the 
cloud 0B7, by Motch et al. (1997) using RASS data 
lend further support to this idea. The survey is com- 
plete to about 0.02 counts"^ and, at this flux level, 
68 sources are detected. Catalogue searches and op- 
tical follow-up observations have shown that the vast 
majority of these sources arc associated with active 
coronae (F-K and M stars). There are 8 NOIDs and 
they do not seem to be correlated with the denser 
regions of the cloud. Again this figure is lower than 
the estimated number of ONSs emitting above 0.02 
counts/s (~ 10 according to Zane et al. 1995). 

A definite confirmation of the paucity of detectable 
ONSs in ROSAT images came from the work of 
Danner (1998a, b). He searched both high-latitude 
molecular clouds and dark clouds in the Galactic 
plane for soft, thermal X-ray sources using RASS 
data. The high-latitude cloud sample used by Dan- 
ner is that originally investigated by Magnani, Blitz 
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& Mundy (1985) and the searched area is ^ 125 deg^, 
substantially larger than that occupied by the clouds 
themselves, ~ 64.5 deg^. The total number of sources 
within the cloud boundaries is 89 and the survey 
should be complete to 0.012 counts"^. Of the 89 
sources 54 are firmly identified with AGNs or stars, 
while for those remaining one, or more, plausible 
counterparts are reported. The claim is that at most 
one source could have escaped the identification pro- 
gram, so the derived upper limit on the projected 
density of ONSs is ~ 0.016 deg~2. This fi gure IS 
smaller than the number expected, ~ 0.15 deg^^ ac- 
cording to Zane et al. (1995), by at least one order 
of magnitude. In his dark cloud sample, which coin- 
cides with that analyzed by Zane et al. (1995), only 
16 sources met the selection criteria and all of them 
have a firm, or at least plausible, optical counterpart, 
in most cases a bright star. The only exception is 
RX J1856. 6-3754, but no other source with similar 
properties was found. Since the dark cloud sample 
covers an area of ~ 1600 deg^ and only one isolated 
NS candidate is detected, the ONS projected density 
is ^ 6 X 10~^ deg~^ at a count rate of 0.05 counts"^. 
A direct scaling of the results in table 5 of Zane et 
al. (1995) to a fixed (i.e. independent of the cloud) 
threshold of 0.05 count s~^ gives 50 detectable sources 
in the same area, nearly two orders of magnitude in 
excess of the observed upper limit. 

7. ACCRETING VERSUS COOLING NEUTRON STARS 

According to the canonical scenario, soon after it 
is born in a type II Supernova event, a NS is very 
hot (T ~ 10^^ K), has large magnetic field {B ^ 10^^ 
G) and short period (P ~ 0.01 s). However, neu- 
trino emission quickly bring the temperature down to 
~ 10^-^ - 10^ K in about 10^ - 10"^ yr (see e.g. Page 
1998). For the next ~ 10^ yr, the star temperature 
stays roughly constant (at ^ 10^ K) and the NS can 
be observed also in the X-ray band 3jS Si weak, soft, 
thermal source. X-ray emission from radio pulsars 
and from the 7-ray pulsar Geminga has been convinc- 
ingly observed (see e.g. Ogelman 1991, 1995; Foster 
et al. 1996 ; Bignami & Caraveo 1996 for reviews; 
see also Becker & Triimper 1998). Many of these 
sources show a two-component spectrum and, while 
the origin of the hard, non-thermal component re- 
mains uncertain, the soft component almost certainly 
originates from the NS surface. Fitting this thermal 
component with a blackbody gives Tf.ff ~ 10^ — 10^ 
K, with lower temperatures associated with older pul- 
sars (see e.g. Pavlov et al. 1995). 

Because of the typical values of T(,f / , L and of the 
thermal spectrum, a middle-aged, radio-quiet, cool- 
ing object appears quite similar to an accreting ONS 
in soft X-rays. This introduces an ambiguity and the 



association of the isolated NS candidates discovered 
so far with either an old or a (relatively) young ob- 
ject remains a matter of lively debate. It is therefore 
extremely urgent to improve our theoretical under- 
standing of these two classes of sources, looking, in 
particular, for spectral signatures which can enable 
us to discriminate between them. 

In principle, the lack of pulsations in the X-ray 
data could be used as an argument against the pres- 
ence of a high magnetic field. This, in turn, may 
suggest that the field has decayed, favoring an inter- 
pretation in terms of an old NS (see section ^) . Un- 
fortunately, this is far from being conclusive since the 
co-alignment between the rotation axis or the line of 
sight and the magnetic axis would produce no pulsa- 
tions even for high B. On the other hand, field decay 
in isolated objects is highly questionable, so the pres- 
ence of pulsations is even less conclusive in estimating 
the age of the NS. 

At least one candidate, RX J1605. 3-^3249 for 
which X-ray measures cover a substantial time in- 
terval, shows a remarkably constant flux on various 
timescales (weeks to years). Motch et al. (1999) ar- 
gued that this can be taken as a hint in favor of the 
cooling picture, because accretion instabilities should 
lead to X-ray variability on a free-free timescale 
which is ~ 1 yr under typical conditions. Also, ISM 
dishomogeneities on the smallest scales may give rise 
to variations in the accretion rate. At a velocity of 
~ 40 kms~^, an ONS would travel a distance ~ 100 
AU in the time separating the first and last ROSAT 
observation of RX J1605.3-h3249 (~ 7.5 yr). 

The (expected) similarity between accreting (low- 
luminosity) and cooling NSs is not surprising. In both 
cases models are based on a geometrically thin, static 
atmosphere in LTE which has roughly the same ther- 
modynamical properties. The only difference is the 
energy input, either the heat released by the imping- 
ing ions or thermal radiation emanating from the NS 
crust. 

The problem of calculating the spectrum emerg- 
ing from cooling NSs has been widely investigated 
by a number of authors, both for different chemical 
compositions and for low and high magnetic fields 
(see e.g. Romani 1987; Miller & Neuhauser 1991; 
Shibanov et al. 1992; Miller 1992; Pavlov et al. 1994, 
1995; Zavlin et al. 1995; Zavlin, Pavlov & Shibanov 
1996; Rajagopal & Romani 1996; Pavlov et al. 1996; 
Rajagopal, Romani & Miller 1997). These studies 
showed that emerging spectra are not too different 
from a blackbody at the star effective temperature. 
For B ^ lO^'' G, spectra from pure H atmospheres 
exhibit a distinctive hard tail that becomes less pro- 
nounced when the magnetic field is ~ lO^'^ — 10^^ G. 
However, as already mentioned in section 3.2, this 
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feature cannot be taken as characteristic of a cool- 
ing object, since similar conclusions were reached for 
the spectrum emitted by low-luminosity, low-field, 
accreting NSs (Zampieri et al. 1995; Zane, Turolla & 
Treves 1999). 

While a cooling atmosphere is likely to be rich in 
heavy elements which results from the supernova ex- 
plosion and the subsequent envelope fallback, the as- 
sumption of a pure hydrogen composition, although 
crude, is plausible for an accretion atmosphere. In 
this case, in fact, incoming protons and spallation 
by energetic particles in the magnetosphere may en- 
rich the NS surface with light elements (mainly H) 
and, owing to the rapid sedimentation, these elements 
should dominate the photospheric layers (Bildsten, 
Salpeter & Wasserman 1992). As detailed calcula- 
tions have shown, the overall emission strongly de- 
pends on the chemical composition. This led to the 
suggestion that the comparison between observed and 
synthetic X-ray spectra may probe the chemistry of 
the NS crust and, indirectly, shed light on the pow- 
ering mechanism. 

Additional information on the nature of these 
sources may come from optical observations. In fact, 
accretion spectra have been found to exhibit a char- 
acteristic excess over the Rayleigh-Jeans tail of the 
best-fitting (planckian) X-ray distribution, which is 
not shared by cooling models with similar parame- 
ters (see Figure 2; Zane, Turolla & Treves 1999). The 
presence of an optical bump is due to the fact that 
in the outermost layers the heating produced by ac- 
cretion cannot be balanced by free-free cooling. The 
temperature then rises until T ~ 10^ - 10^ K, when 
Compton cooling becomes efficient. The presence of 
this hot "corona" does not influence the propagation 
of X-ray photons, which come from the inner lay- 
ers, but makes the optical photosphere slightly hotter 
than in the cooling case. 

The observation of an optical excess has been re- 
ported in the spectrum of RX J18563. 5-3754 (see sec- 
tion |6.2| ). Although the optical identification of RX 
J0720. 4-3125 still lacks a definite confirmation, it is 
interesting to note that the counterpart proposed by 
Kulkarni & van Kerkwijk (1998) also shows a similar 
excess. 

The multiwavelength spectrum of RX J18563. 5-375 
has been the subject of various investigations. Mod- 
els of cooling atmospheres based on different chem- 
ical compositions fail to predict the correct optical 
fluxes (Pavlov et al. 1996), while models with two 
blackbody components or with a surface temperature 
variation may fit the /3000 flux. In particular, the fact 
that a simple model with T oc (latitude) '^■^^ and with 
magnetic and rotational axes co-aligned can match 
the observed spectral energy distribution (SED), fa- 



vored the interpretation of RX J18563. 5-3754 as a 
cooling NS (Walter & Matthews 1997). Recent spec- 
tra from non-magnetic atmospheres with Fe or Si-ash 
compositions (Walter & An 1998) may also provide 
a fit of both the X-ray and optical data, although 
they still need a definite confirmation. The magne- 
tized model by Zane, Turolla & Treves (1999) indeed 
predicts an optical excess of about the right order in- 
dicating that, given the considerable latitude of the 
unknown parameters, the full SED may be consistent 
with the picture of an accreting NS. However, no de- 
tailed fit for RXJ 185635-3754 has been attempted 
yet. It should also be noted that accretion spectra, 
irrespective of details on their shape, are intrinsically 
harder because of the reduced emitting area (see eq. 
p5[|). Assuming / ^ 0.01, which is conservative if the 
field is high and equation ([l^) is taken literally, the lu- 
minosity corresponding to a peak energy of ~ 100 eV 
IS 1028-1029 ergs ^s, more than one order of magni- 
tude lower than what follows from observations plac- 
ing the source at ~ 100 pc. 

There has been a preliminary report of a proper 
motion measurement for RXJ 185635-3754 by the 
MPE group (Triimper, private communication; Fred 
Walter's web page). If confirmed, the value of ~ 
0.34 arcsec/yr would imply a projected velocity of 
~ 160 (d/lOOpc) kms^^. Such a high velocity is def- 
initely incompatible with the observed fiux, if the 
source is powered by accretion and d ~ 100 pc. Nev- 
ertheless, the accretion option could still be consis- 
tent with the data placing the source at a closer 
distance. For (i ~ 30 pc, a velocity ~ 50 kms~^ 
would give both the correct values of the proper mo- 
tion and the luminosity. Besides, a distance of the 
same order was derived by Campana, Mereghetti &: 
Sidoh (1997) from the fit of ROSAT data with the 
unmagnetized synthetic spectrum of Zampieri et al. 
(1995). Two major objections still remain. First, the 
observed column density is too large to be consistent 
with d ~ 30 pc, taking into account that the local 
ISM is underdense. Second, the probability to find 
a NS so close and with such a low velocity must be 
extremely small, also in the light of the recent result 
by Popov et al. (1999). Here we just mention that, 
at variance with cooling atmospheres, all accretion 
spectral models were computed so far under the as- 
sumption of complete H ionization and do not include 
the contribution of the atmospheric neutral atoms to 
the column density. A rough estimate shows that 
a small fraction of neutrals (^ 10%) is anyway ex- 
pected and this may produce the required absorption 
"in loco". As far as the second point is concerned, it 
is worth noting that the chance of finding a slow ONS 
at 30 pc is not much slimmer than that of seeing a 
fast, cooling NS at 100 pc. Taking a total number of 
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NSs ~ 10^, the average ONS separation is ~ 15 pc, 
so the probabihty of finding it at 30 pc is essentially 
1. According to Popov et al. (1999) the fraction of 
low-velocity NSs is less than 1%. This gives a total 
probability of ^ 1%. Cooling NSs are short-lived, 
tcooi ~ 10^ — 10® yr (Page 1998) and their average 
separation is larger, ~ 270 pc. Now the probability 
of finding a cooling NS with V ~ 240 kms~^ is es- 
sentially unity, but the chance of having it inside 100 
pc is therefore ~ (100/270)^ ~ 0.05. The conclusion 
is that both events are quite rare. 

The problem of distinguishing between accreting 
and cooling isolated NSs on the basis of a statisti- 
cal analysis has been recently faced by Neuhaiiser & 
Triimper (1999). They have shown that the logA^- 
log S curve for the observed candidates is closer to 
that of coolers, at least if ONSs are taken as a slow 
population as in TC and BM, and suggested that 
higher mean velocities may yield better agreement 
with observed upper limits. This particular point will 



be addressed in section 8.2. 



8. WHY ARE SO FEW NSS DETECTED? 

As discussed in section 7, the first analyses of 
ROSAT survey data are mainly negative as far as 
ONSs are concerned, leading to a problem of paucity 
of detections. The theoretical predictions of the "pre- 
ROSAT" epoch appeared to be quite optimistic, forc- 
ing theoreticians to reconsider critically their starting 
assumptions. While the ONS spectral properties are 
rather well established, earlier statistical estimates 
rely on a number of crucial assumptions concerning 
the velocity distribution of pulsars at birth and the 
long term evolution of the magnetic field. Previous 
studies (see section 5.2) were based on the hypoth- 



esis that nearly all ONSs have a residual magnetic 
field ~ 10^ G and are presently in the accretor phase. 
However, accretion is just one of the different phases 
that a NS can experience. Normally born in the 
ejector stage, the NS enters later the propeller, and 
the accretor phase if the conditions are favorable. A 
time-dependent magnetic field may therefore change 
the duration of these phases and alter the statistics 
of visible ONSs. The evolution of the i3-field in iso- 
lated NSs is still a very controversial issue. Little 
evidence comes from the pulsar statistics, and one 
should rely only on theoretical models (see e.g. Srini- 
vasan 1997). Theoretical results are far from being 
univocal and predict either exponential/power-law 
field decay (Ostriker & Gunn 1969; Sang & Chan- 
mugam 1987; Urpin, Chanmugam & Sang 1994; Miri 
1996; Urpin, & Muslimov 1992; Urpin & Konenkov 
1997) or little or no decay at all within the age of the 
Galaxy (Romani 1990; Srinivasan et al. 1990; Gol- 
dreich &: Reisenegger 1992; see also Lamb 1992 for 



a review). Statistical analyses based on observations 
of isolated radio pulsars are consistent with no field 
decay over the pulsar lifetime (Narayan & Ostriker 
1990; Sang & Chanmugam 1990; Bhattacharya et al. 
1992). However, this does not preclude the possibil- 
ity of field decay over longer timescales. Different 
approaches to pulsar statistics led, independently, to 
the conclusion that, if the magnetic field decays, then 
it probably does so over a timescale ~ 100 Myr (Hart- 
mann et al. 1997). 

A step forward in the statistical analysis is there- 
fore to assume different laws for the time evolution of 
the field and to combine the dynamical evolution of 
the simulated population with its magneto-rotational 
evolution. This has been computed for the first time 
by Popov et al. (1999) and is discussed in section 9.2. 
In section 9.1, we briefiy summarize the role played 
by the magnetic field in affecting the NS magneto- 
rotational evolution. 

It should also be mentioned that the Bondi-Hoyle 
accretion rate (eq. fl^) may be just an upper limit. 
As discussed by Blaes, Warren & Madau (1995) the 
EUV-X radiation released in the accretion process 
can heat the surrounding medium, producing an in- 
crease of the sound speed which, in turn, tends to 
decrease M. This effect, known as preheating, may 
severely reduce the accretion rate in ONSs. 

8.1. The role of the magnetic field 

The permanence of a NS in the ejector (or pro- 
peller) phase depends on the magnetic field strength 
and is influenced by its time evolution. Livio, Xu & 
Frank (1998) and Colpi et al. (1998) were the first to 
recognize that a NS can straggle in the ejector or pro- 
peller state for a time longer than the present cosmic 
time mainly as a consequence of magnetic field decay. 
This may cause the paucity of accreting ONSs. On a 
theoretical basis, the decay occurs either because of 
ohmic dissipation of crustal currents sustaining the 
field (Sang & Chanmugam 1987; Urpin & Konenkov 
1997), or because of the migration of the proton flux- 
oids from the core to the crust due to continuous 
spindown (Srinivasan et al. 1990). 

The idea that a NS may never become an accre- 
tor is rather simple and can be illustrated by just 
estimating the duration of the ejector stage te and 
introducing a simple exponential decay law for the 
field. Magnetic dipole braking in this phase keeps 
the neutron star spinning down at a rate 



dipole 



Bjt) 
1012G 



p-^s yr"^ 



(20) 



until Pcrit is reached (see eq. [6] ) . The critical period 
is itself a function of time [Pcrit oc i?^/^) and de- 
creases as B{t) decays. Two competing effects thus 
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come into play: (i) the spindown rate slows down be- 
cause of the weakening of the field, causing the NS 
to stay longer in the ejector state; (ii) Pcrit decreases 
with time and this acts in the opposite way, making 
te shorter. If the field decays exponentially on a scale 
Td from its initial value Bq down to a minimum value 
i?b (the presence of a relic field is mainly suggested by 
the stability of the field in millisecond pulsars), the 
duration of the ejector stage can be computed analyt- 
ically. Following closely Popov & Prokhorov (1999), 
the ejector time scale is 



TE 



n + 




TE < Tb 



TE > n, 



(21) 

where Tf, = ln(i?o/-S;,) is the time at which the NS 
reaches the minimum field and 



TO ~ An-^/\io{Bo/W^^G)-^ Gyr 



(22) 



is the duration of the ejector phase for a constant 
field; note that tq depends on v and n. Figure (^) 
shows the loci where the time spent in the ejector 
phase equals the age of the Galaxy, te = 10 Gyr, 
as a function of the minimum field and of the de- 
cay timescale. Different curves refer to initial fields 
of 5 X 10l^ 10^^ 2 X 10^2 G and u = 40 kms'^, 
n = 1 cm~^ were assumed. The "forbidden" region, 
where r^; > 10 Gyr, lies below each curve. As is ap- 
parent from figure ^ the parameter range for which 
the NS never leaves the ejector stage is significant 
even for the rather low value of the star velocity we 
adopted. It is interesting to note that unhalted decay 
(i.e., Bh — > 0) would drive all NSs into the accretor 
phase, since the Alfven radius becomes smaller than 
the star radius as the field decays. 

Colpi et al. (1998) considered a model in which the 
spin evolution causes the core field to migrate to the 
crust where dissipation processes drive ohmic field de- 
cay. They found that for a characteristic time scale 
of ~ 10* yr (a value consistent with the limits on the 
decay time scale inferred from pulsar statistics) the 
number of ONSs that are at present in the accretion 
stage can be reduced by a factor of 5 over previous 
figures. Most of the (low velocity) ONSs would now 
be in the ejector or the propeller stages. 

8.2. Are ONSs Fast ? 

Pulsars have large kick velocities. Thus, if NSs 
at birth acquire large mean speeds their distribution 
may be devoided of low velocity stars. This may ex- 
plain the paucity of observed ONSs. In the light of 



new data on pulsars (see § |]^) a re-analysis of the 
statistical properties of NSs is required. 

Popov et al. (1999) studied the effect of an increas- 
ing initial mean velocity on the visibility of ONSs. 
In this investigation, the NS spin-down induced by 
dipole losses and their interaction with the ambient 
medium (described using a realistic map of the ISM 
in the Galaxy) has been followed, together with their 
dynamical evolution in the Galactic potential, in un- 
precedented detail. Not unexpectedly, the ejectors 
were found to outnumber the accretors (and the pro- 
pellers): they comprise more than 95% of all NSs at 
(V) ~ 150 kms-^ The visible ONSs fah off rapidly 
when (V) is above 200 kms~^, owing to the lack of 
low velocity stars with increasing mean speed. 

Using the present upper limits on the number of 
detected ONSs Popov et al. (1999) were able to 
constrain the pulsar velocity distribution indepen- 
dently of radio observations. Taking a (generous) 
upper limit of 10 accreting sources within 140 pc 
from the Sun, these authors inferred a lower bound 
on the mean kick velocity (V) of 200-300 kms~^, 
corresponding to a dispersion ay ~ 125-190 kms~^, 
in agreement with the findings of Cordes & Cher- 
noff (1998). In addition, these results have been able 
to constrain the fraction of pulsars in the low veloc- 
ity tail which could have escaped pulsar statistics, to 
< 0.05%. 

9. CONCLUSIONS 

In the last decade old neutron stars accreting from 
the interstellar medium have been actively sought, 
mainly in X-ray ROSAT images. Although the quest 
produced half a dozen strong candidates, none of 
these can be unambigously claimed as an accreting, 
isolated NS, so that the very existence of this class of 
sources is not observationally proven as yet. 

The expected number of observed ONSs depend on 
the physical properties of their progenitors, young 
NSs in the pulsar phase. The key quantities are 
their velocity distribution and the law of magnetic 
field decay, both subjects beeing frontline problems 
in pulsars physics. Pulsars and ONSs would be the 
A and 17 of the evolution, which brings the neutron 
star through the intermediate phases of ejector and 
propeller. 

The physics of propellers is poorly known and, es- 
pecially in the case of very low accretion rates, it 
is suspected that emission could be non-stationary. 
If flaring episodes do occur, they may be of critical 
importance in recognizing a population, which is oth- 
erwise inaccessible because of its intrinsic low lumi- 
nosity. 

Substantial improvements of our current under- 
standing of the isolated NS candidates discovered 
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so far are expected from next generation X-ray 
satellites. Advanced spectroscopy with CHANDRA, 
XMM, and ASTRO E will provide high-resolution 
X-ray data which can reveal many properties of the 
NS atmosphere, like composition, gravitational red- 
shift, temperature. The comparison of high-quality 
spectra with models might indeed prove decisive to 
assess the nature of these sources. Optical observa- 
tions with the largest telescopes, already operational, 
like KECK and VLT, would prove of the same im- 
portance in solving the riddle of isolated NS candi- 
dates. Without the identification of optical counter- 
parts they will remain X-ray bright NOIDs. High- 
quality spectroscopy of the optical counterparts is 
needed to assess better the (possible) deviations from 
the Raylcigh-Jeans tail of the X-ray spectrum, like 
those reported in two candidates so far. 

Hopefully the next decade will bring us strong new 
candidates for isolated NSs steadily accreting from 
the ISM, and the discovery of the previous unsteady 
phase. 
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Fig. 2. — Synthetic spectra of pure H accreting (upper curves) and cooling (lower curve) NSs. a): B = G, L — 4.3 x 10'^"^ crg/s; 
h): B = 10^^ G, L = 3.7 x 10^^ erg/s; c): B = 10^^ G, L = 0.2 x 10^^ erg/s. The blackbody spectra which best-fit the X-ray 
range close to the peack of each spectrum are also shown (dashed curves). The two arrows mark the optical band (3000-7500 A). 
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